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Abstract 

We construct charged anisotropic AdS domain walls as solutions of a consistent truncation 
of type IIB string theory. These are a one-parameter family of solutions that flow to an AdS 
fixed point in the IR, exhibiting emergent conformal invariance and quantum criticality. 
They represent the zero-temperature limit of the holographic superfluids at finite superfluid 
velocity constructed in larXiv: 1 010.5777, We show that these domain walls exist only for 
velocities less than a critical value, agreeing in detail with a conjecture made there. We also 
comment about the IR limits of flows with velocities higher than this critical value, and point 
out an intriguing similarity between the phase diagrams of holographic superfluid flows and 
those of ordinary superconductors with imbalanced chemical potential. 



1 Introduction 



Quantum critical points are expected to be of significance in understanding the ground 
states of high-T c superconductors. Holographic constructional of such quantum critical 
(hence zero-temperature) superconductors give rise to domain wall solutions, which cap- 
ture the holographic RG flow from a symmetric state in the UV to a symmetry-breaking 
vacuum in the IR pp. 

Examples of such solutions were presented in J3J, 0]. They arise in consistent truncations 
of type IIB and M-theory containing an abelian gauge field and a charged scalar with a 
symmetry-breaking potential. To trigger the RG flow, the UV theory was deformed by 
a chemical potential fi. That is, the U(l) symmetry is not broken explicitly, rather the 
theory is deformed by means of an uncharged operator /iJ° (i.e., Aj^J 11 with only the time 



component turned onil). This induces a non- vanishing VEV for the charged scalar, hence 
breaking the U(l) symmetry spontaneously. The domain walls interpolate between two AdS 
minima, in the UV and IR, and represent the ground states of holographic type IIB and 
M-theory superconductors constructed in [5] and [6]. 

A natural generalization of this idea is to break the isotropy by turning on a spatial 
component of the boundary vector potential. The corresponding solutions can be interpreted 
as the zero temperature limits of holographic superfluid flows (i.e. a superfluid with finite 
superfluid velocity £) discussed first in [7J |8] and generalized in various ways in [91 [101 fTTj . 
In [T2], fully backreacted holographic superfluid flow solutions were constructed at finite 
temperature within the same type IIB consistent truncation of [3J. The advantage of a 
backreacted solution was that one could systematically lower the temperature in a fully 
controlled way. Using this approach, strong indications were obtained in [12] that in the zero- 
temperature limit the IR AdS space found by [3J is actually robust even when a superfluid 
velocity is turned on, at least up to a critical value £ c . Various pieces of evidence were 
presented that there should exist an anisotropic domain wall solution interpolating between 
two AdS spaces, for velocities in the range < £ < £ c . On the other hand, the holographic 
superfluid flows constructed in [12] exist in the bigger range < £ < £*, with > £ c being 
the velocity above which superfluidity is destroyed and the system is always in the normal 
phase. The low temperature behavior of the solutions in the range £ c < £ < indicated 
that a different non AdS-like IR phase emerges at zero temperature. In this range the 
perturbation induced by £ turns out to be too strong to be washed out by the RG-flow, so 
that the anisotropy still survives at small radii, suggesting the non-existence of any obvious 
quantum critical point. 

1 Holographic superconductors were first constructed in [2]. 

2 The superscript B denotes the fact that these are boundary quantities. 




1 



In the present paper we put the speculations of [12] on a concrete footing by explicitly 
constructing gravitational solutions describing the ground state of holographic superfluid 
flows. We find that anisotropic AdS-to-AdS domain wall solutions indeed exist only in the 
expected range, < £ < £ c , and make various quantitative checks that these are the zero- 
temperature limits of the superfluid flows constructed in [12]. On the contrary, for higher 
velocities we do not see the emergence of an AdS-like geometry in the IR. It is an interesting 
question to ask whether there is a meaningful IR limit in this regime, and if so, whether full 
solutions can be constructed explicitly. While we did not succeed in finding any meaningful 
geometry emerging in the IR for £ c we offer a few preliminary comments in the 

concluding section. 

T A 



Quantum 
criticality 




Figure 1: The (qualitative) phase diagram of holographic superfluids. At zero temperature, a 
quantum critical point is found for velocities below a critical value, £ c . Above £ c the system 
enters a more anisotropic phase where the deformation induced by £ affects the RG-flow 
strongly, and brings the system away from any obvious AdS-like IR fixed point. 



Intriguingly, the phase diagram which emerges from our analysis, Figure [H is reminis- 
cent of those expected for real-life superconductors when one induces (e.g., by an external 
magnetic field) an imbalance 5fi in the chemical potential between the two populations 
of fermionic particles that form Cooper pairs^l For high enough S/j,, still lower than the 
critical value above which superconductivity is completely destroyed, the system enjoys an 
anisotropic phase at T = known as FFLO phase. This phase is characterized by a sponta- 
neous breaking of translational invariance produced by the imbalance <5/i, and is separated 
from the ordinary superconducting phase by a Chandrasekhar-Clogston (CC) bound. In this 
regime, the system finds it energetically favorable to be in a superconducting anisotropic con- 
figuration, where Cooper pairs have a non zero net momentum p3]. While this has been 
proven for BCS superconductors [T5], it is argued to be a generic phenomenon also holding 

3 We thank A. Cotrone for bringing our attention to these systems. 
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e.g. in high-T c superconductors. Although the physical setups are different, the response of 
this system to S/j, is very similar to the response of our holographic setups to a non- vanishing 
superfluid velocity £. It would be interesting to investigate this analogy further, and more 
generally, use holography to study superconductors with an imbalanced chemical potential, 
including the prediction for a CC bound. 

The rest of this paper is organized as follows. In section 2 we present the effective 
five-dimensional bulk action, our ansatz for the various fields, and discuss the IR and UV 
boundary conditions our solutions should satisfy. Section 3 contains our one-parameter 
family of solutions and a discussion of their properties. Finally in section 4, we comment 
on our results and discuss the regime for which AdS domain wall solutions do not seem to 
exist. 



2 The model: construction and strategy 

The theory we will use is the consistent truncation of low energy type IIB string theory 
considered in [5]. The action takes the form 



Sub = / d 5 xy/^g 



L 2 1 /2/\ 3 

R - —F ah F ab + - (- e abcde F ab F cd A e + 



1 I -/-\2 , -1.2 j./a a o A \2 ^ 2 ( if) 



- I (d a ipy + sinh 2 ^(d a 6 - 2A a ) 2 - — cosh 2 ( £ J (5 - cosh vo 



• (2-1) 



We have set I&ttG = 1 and the conventions are such that e 01234 = 1 / y/—g, and we have 
written the charged (complex) scalar by splitting the phase and the modulus in the form 
ijje . The Abelian gauge field A is dual to an i?-symmetry in the boundary field theory [5] 
and the scalar field has i?-charge R = 2. The general equations of motion that follow from 
this action (working in the gauge 6 = 0) can be found in [12] . 

In [T2] , strong evidence was presented for the existence of an anisotropic AdS domain wall 
that was a solution of the above Lagrangian. Our aim here is to construct it explicitly. To 
do this, we will write down an ansatz for the metric, gauge field and scalar with only radial 
dependence, so that the resulting equations of motion are ordinary differential equations. 
We will take the same ansatz as was used in [T2]. This is equivalent to an anisotropic 
generalization of the structure considered by Gubser et al. in [3]. Specifically, we take 

ds 2 = - r -^f-dt 2 + ^P-dr 2 - 2C(r) r -dtdx + ^-B(r)dx 2 + "-dy 2 + T -dz 2 (2.2) 
L A ij(r) L z L z L A L z 

for the metric, and 

A = A t (r) dt + A x (r) dx , ip = i/j(r) , (2.3) 
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for the gauge and the scalar fields. 

With this ansatz, the full set of equations of motion can be massaged to the form of 
two first order and five second order differential equations for seven functions. Since the 
structure of the ansatz is the same, the equations of motion we get are also the same as in 
[T2"] . As discussed there, these equations exhibit four scaling symmetries and we quote them 



here for convenience 

£->-t/a, / ->■ a 2 / , h^-ah, C ->■ aC , A t -ta,A t , (2.4) 

x^x/b, B^b 2 B, C^bC, A x ^bA x , (2.5) 

(r,t,x,y,z,L) -> a(r,t,x,y, z, L) , (A t , A x ) (A t , A x )/a , (2.6) 

r^/3r, (t,x,y,z) (t,x,y,z)/P , (A, A.) -> /3(A t) A,) . (2.7) 



We can forget the third scaling symmetry by setting the length scale L — 1. The remaining 
scalings will be very useful in what follows. 

So far our setup is identical to that in [12]. The crucial difference is that the IR boundary 
conditions we will impose now are those of an IR AdS space, not a black hole horizon. In 
our coordinates, IR corresponds to r — > and UV to r — > oo. Our aim now is to construct 
a domain wall solution that interpolates between two AdS vacua, corresponding to the flow 
between the symmetry-preserving (UV) and symmetry-breaking (IR) vacua of the scalar 
potential. To compute the IR expansion, we need to know the IR AdS background. The IR 
AdS scale is fixed by the value of the scalar potential at its symmetry-breaking minimum in 
the action above. The minimum is at 

= ArcCosh 2 , (2.8) 
and the effective AdS scale there can be computed from the action to be 

2 3/2 

L IR = —L. (2.9) 

Unlike the finite temperature case, the correct IR boundary expansions of the fields are now 
no longer of a simple power series form in r. The easiest way to see this is to note that the 
scalar field equation in the IR AdS background allows the solution 

ip = ArcCosh 2 + ^i,o r a , (2.10) 

with ipi t0 an integration constant. The reason for choosing this specific form for the subscripts 
on this coefficient will become clear momentarily. Here, a can be computed^) by solving the 
scalar field equation in the IR AdS, and the result is 

a = 2^3- 2 ~ 1.46. (2.11) 
4 This value was computed also in [3], but the value reported there has a typo. 
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Since we are looking for a non-trivial solution that starts off at the symmetry-breaking 
minimum for small r, we should allow for such irrational powers in the IR expansion when 
trying to extrapolate from the IR to UV. We will try out an expansion^] for the fields by 
writing them in powers of the two variables x and y defined as 

x = r a , y = r, (2.12) 

and then constrain the coefficients by demanding that they satisfy the equations of motion. 
This fixes the form of the fields in the IR (up to 0(r 2 )) to be 



f(r) = ^0,0 + A -^f^ l> 2 + - , Hr) = h , + h 0<2 y 2 + ... , (2.13) 

B{r) = B , - y 2 + ... , C(r) = C o , + 4A ^f^ + ... , ( 2 .14) 

A t (r) = A m) y 2 + - , A x (r) = A x{0 , 2) y 2 + ... , (2.15) 

ip(r) = ArcCosh 2 + Vi,o x + ^0,2 y 2 + ••• , (2.16) 

where not all of the coefficients are independent. In particular 

_ 2(-27/^ A 2 x{ ^ 2) + 64Cff [0 A 2 { ^ 2) + 48/^ C 0fi A t{0j2) A x{0t2) + 96/^ B 0j0 A 2 { ^ 2) ) 
0,2 ~ 27(8C 2 )0 /io,o + 95 ,o^,o) ' 

(2.17) 

lQ(9h 2 A 2 x(0)2) - 16C ,0 A(0,2) Ar(0,2) - 8^0,0 ^(0,2)) n oN 

Wno = = . (Z.lo) 

3V3(8Cg,o + 9Bo,o^ i0 ) 

The notation for the coefficients should be clear by now: in the IR expansion for a field 
<f), the coefficient of x a y b is denoted as (f) a ^. All higher order coefficients (which include 
powers of x, powers of y and mixed powers) are fixed iteratively in terms of the lower order 
coefficients due to the equations of motion. We will need the IR expansion to start the 
numerical integration from the IR, but the explicit forms of the higher coefficients are not 
very enlightening, so we will not report them here. Note, finally, that there are only six 
independent quantities in the IR, just as there were at the horizon in [12J. In the present 
case, these independent quantities are 

{^o,o, B 0fi , C 0)0 , At(p,2), A x (p >2 ) , ^1,0} • (2.19) 

Our strategy in constructing the solutions will be to pick numerical values for these IR 
data, and integrate the equations of motion numerically all the way to some very large r, 
corresponding to the UV. At the UV boundary, we can again solve the equations of motion 



3 We thank S. Gubser and S. Pufu for helpful discussions on this. 
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straightforwardly and look for a consistent expansion that reproduces the curves arising 
from this integration. The relevant series expansions can be found in eqns. (3.18)-(3.23) 
and Appendix A of [12] and we will use the notations there in our UV discussion. To get 
asymptotically AdS boundary conditions in the UV, we need to set B = 1 = h , and 
Co = = ipi. The former two conditions can be accomplished via the rescalings (12.41 
12.51) . while a shooting technique is required for the latter two. This gives rise to the eight 
independent boundary quantities [12J: 

(U , B 4 , C 4 , A t;0 , A t! 2 , A Xi0 , A Xj2 , Vs) • (2.20) 

We will also use the rescaling (12 .7p to set the leading piece of A t at the boundary (namely 
A t) o = fi) to unitj@. Note that since this rescaling involves r, we will need to appropriately 
rescale the coefficients of the IR expansion as well, when launching the integration^. Once 
we fix the chemical potential to one, for any given superfluid velocity £ (before the rescaling 
of the solution that sets A t to unity, £ is given by A Xt0 /A t fi) the number of independent pa- 
rameters at the boundary is six, which is the same as the number of horizon data. Therefore 
we expect to find at most discretely many domain wall solutions for any given superfluid 
velocity. As discussed in [3] we will look for the solution where the radial profile of the scalar 
field has the least number of nodes. 



3 Results 

After the shooting and the rescaling discussed in the previous section, we have the domain 
wall solution we were seeking. We present the plots of the various functions for a selected 
value of the velocity in Figure 2. We have checked that the plots for other velocities are 
qualitatively similar. We find that solutions exist for velocities £ ranging from to £ c 0.374, 
which is consistent with the expectation of [T2]: there we found evidence that an AdS IR 
fixed point exists for low enough velocities but not for £ > 0.40, even though the finite- 
temperature superfluid phase exists to a higher value ~ 0.5. It was also found there that 
the condensate, as defined by 

{ODw) = (3.1) 

tends to a fixed value as the temperature is lowered, for small enough velocities. (In writing 
the above expression for the condensate, we use the fact that A t $ has been rescaled to one.) 

6 In [5] this scaling was used to set the horizon datum ^t(o,2) = 1- We prefer instead to set At,o = 1 at 
the boundary: this corresponds to working in a fixed chemical potential ensemble in the gauge theory. 

In the finite temperature case of [12] . this scaling shifts the horizon radius and therefore effectively 
introduces a new parameter, the temperature. 
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Figure 2: Plots of the various functions as a function of the radial coordinate for £ = 0.33. 
A logarithmic radial coordinate has been chosen to illustrate that the solution is a domain 
wall. For other values of the superfluid velocity in the range < £ < £ c the nature of the 
curves is similar. 



We can compute the same quantity using our domain wall solutions, and we find that for all 
the velocities for which it exists, the value of this condensate is a constant and is precisely 
equal to the one found as the zero-temperature limit of our finite temperature solutions. We 
can also perform another check of our solutions, by considering the isotropic (i.e., £ = 0) 
limit. In that case, it was found in [3J that the object (O)gpr = ( 2 /t/3) 3 ~ 0-322. We have 
checked that our solutions reproduce this value in the isotropic limit. 

The solutions we constructed were obtained via a shooting method. So one might worry 
that the absence of solutions above £ ~ 0.374 is an artifact of the numerics. It is difficult to 
disprove this conclusively, but we can look at the profiles of some quantities as the velocity 
is changed and see what happens to them as the critical velocity £ c is approached. One 
such useful quantity is the refractive index between the UV and the IR. We can define this 
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quantity as the ratio of the propagation velocities of light in the UV and IR 



n 




(3.2) 



where / is contained in the time component of the metric ( 12. 2p . We plot its behavior as 
a function of the superfluid velocity in Figure 3. It is clear from the figure that at £ c it 




Figure 3: Plot of the refractive index as a function of the superfluid velocity. The refractive 
index diverges around £ c ~ 0.374, where the IR AdS space collapses and the domain walls we 
have constructed cease to exist. In the isotropic limit the numerical value of the refractive 
index is n w 2.674 which is identical to what is reported in pQ. 



diverges, strongly suggesting that domain wall solutions of this kind do not exist beyond 
£ c . It turns out that the AdS space in the IR degenerates as £ — > £ c because /i ,o m the IR 
expansion ( I2.13P tends to zero. 

4 Interpretation and discussion 

Let us summarize our results and discuss them in light of the Criticality Pairing Conjec- 
ture (CPC) proposed in [3]. We can state the CPC as follows. Consider a theory with a U (1) 
symmetry that, once deformed by some appropriate operator, flows to an IR fixed point that 
breaks the U(l). Then the claim is that the theory or its deformations by operators that 
do not break the U(l) explicitly has a finite density, zero temperature state whose infrared 
behavior is controlled by the same IR fixed point. 
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The deformation considered in [3] is given by fiJ°. This operator doesn't break the 
£7(1) symmetry, and hence satisfies the premise of the CPC. The explicit construction of the 
domain wall in [3] demonstrates that the deformation by this operator results in a flow to 
AdS in the IR. In this paper, we considered a more general deformation, fiJ° +A Xi qJ x , which 
breaks isotropy but which again does not break the U(l) symmetry. We find that when this 
more general deformation is added, the system is controlled by the same AdS IR fixed point, 
in agreement with the CPC. However, this holds only as long as the source term, A x< q, is 
small enough: for high enough velocity we seem to find a counter-example to the CPC. 

A violation of the CPC was presented in [TB] where for small enough values of the charge 
of the scalar the IR fixed point was destabilized. In our case, the parameter is not directly 
visible in the supergravity Lagrangian, and therefore the destabilization that happens at 
large velocities is of a different kind. Note that since we are working at zero temperature 
and the theory is conformal, all non-zero values of the chemical potential are equivalent. 
But when we turn on A x ^ as well, we have a tunable parameter, namely the superfluid 
velocity £ = A x $j 'fx. Once this becomes too large, the anisotropy becomes too strong to be 
washed out in the IR. This is similar to what happens in ordinary superconductors: when the 
chemical potential imbalance becomes large, they enter an anisotropic FFLO phase [HI [To] . 

In fact, one of the questions we have left unanswered in this paper is that about the 
possibility of constructing a useful zero-temperature limit for the superfluid flows in the 
range £ c < £ < £*. As already mentioned, in [IB] it was found that the CPC can be violated 
if the scalar charge is tuned to be small, and in such case a Lifshitz geometry emerges in the 
infrared. This is not unreasonable because spatial isotropy is not broken in that case - only 
a charge density is turned on, and no current. In our case, spatial isotropy is broken^ and 
our preliminary investigations suggest that the system is not controlled by a simple Lifshitz 
(or generalized Lifshitz [T7]) geometry. This can be understood as follows. A geometry of 
the form 

ds 2 = -r Zt dt 2 + r Zr dr 2 + r Zx dx 2 + r Zy dy 2 + r Zz dz 2 (4.1) 
has curvature scalars R and R a b c dR abcd that behave as 

R , RabcdR abcd ~ • (4-2) 

(Generalized) Lifshitz corresponds to the case z r = —2 [17], and then both quantities are 
constants. But numerically, as we lower the temperature of the finite temperature solutions, 
we find that there is no region close to the horizon where both R and R a bcdR abcd stabilize 
simultaneously. This could of course be a limitation of the numerics, but together with 
the fact that this is happening at high velocities is a suggestion that if an IR geometry 



? As we have seen, this is not a sufficient condition for the IR geometry to be anisotropic, but it is a hint. 
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exists, its anisotropy should take a different form. Allowing for z r 7^ —2 might be another 
possibility^, so is the possibility of allowing non-trivial scalar and gauge field profiles in the 
IR. Note, however, that just allowing non-trivial matter fields will not bypass the constraints 
on the curvatures (14. 2p . The bottom line is that it would be interesting to construct such an 
anisotropic AdS-to-nonAdS domain wall, if it exists. Our preliminary attempts to look for 
one have been inconclusive. As reported in [12] , at finite but small temperature the curvature 
scalars begin to grow quite quickly as we approach the horizon, so it is not obvious that there 
is a sensible way in which one can assign an IR geometry to this flow. It is possible that the 
natural interpretation is that for velocities greater than £ c the superfluid exhibits runaway 
in the IR and not a fixed point. Answering this question might also let one understand to 
what extent the analogy between holographic superfluids at finite superfluid velocity and 
superconductors with imbalanced chemical potential, can actually be pushed. 
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